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1.  Introduction 

A  great  deal  of  work  has  been  and  is  being  dene,  many  papers  have  been 
published,  and  many  compvitur  programs  are  available  concerned  with  calculating 
pressure  distributions  over  lifting  wings  of  finite  aspect  ratio.  Nevertheless,  the 
state  of  the  art  is  not  satisfactory.  New  approaches  continue  to  be  proposed,  not 
only  because  of  various  complications  in  the  wing  model,  but  largely  due  to  inherent 
difficulties  which  arose  from  singularities  in  the  mathematical  formulation  of  the 
problem.  Consider  here  only  the  basic  engineering  model,  the  thin  lifting  surface 
in  steady  linear  subsonic  (hence,  by  the  Prandtl-Glauert  transformation,  incom¬ 
pressible)  flow.  No  exact  solution  for  any  problem  of  this  typo  lias  been  available. 
Available  numerical  methods  are  approximative  collocation  methods.  The  agreement 
between  calculated  result  and  experiment  is  often  satisfactx  y;  on  the  other  hand, 
there  arc  exceptions.  Then,  in  the  absence  of  an  exact  sedition  for  the  necessarily 
idealized  model  of  the  analysis,  there  is  no  criterion  by  wl  *ch  to  determine  if  the 
fault  lies  with  either  the  idealizations  in  the  model  or  the  inaccuracies  of  the  colloca¬ 
tion  method,  The  temptation  is  to  blame  the  model,  but  div  ergencies  continue  to 
occur  in,  and  betv/oen,  collocation  analyses  for  a  given  mo-del.  Hence  efforts  to 
improve  collocation  analysis  continue. 

The  present  paper  docs  not  deal  with  collocation  analyses.  Its  purpose  is  to 
make  exact  f-olutions  readily  available  - -rot  exact  solutions  for  all  possible  engineering 
models,  of  course,  but  exact  solution  for  the  basic  model  in  the  shape  of  a  circular 
wing.  These  solutions  describe  the  pressure  singularity  th?.t  arises  on  a  parabolic 
tip  of  a  wing  of  any  aspect  ratio  (one  of  the  hitherto  unresolved  fundamental  problems 
of  potential  flow  theory)  but  they  are  useful  beyond  this  somewhat  academic  point. 

They  cv  be  used  to  chock  the  reliability  of  a  given  collocation  method.  They  can  also 
provide,  sometimes  by  means  of  a  very  brief  and  straightforward  calculation,  reliable 
insights  into  the  mechanism  of  lifting  flow  which  would  bo  very  difficult  and  cumber  - 
some  to  derive  convincingly  with  a  collocation  analysis.  To  an  extent,  they  can  be 


2 


k 

generalized  to  other  wing  planforms.  Furthermore,  they  allow  one  to  calculate 
the  initial  vortex  trail  exactly. 

Our  exact  solution  is  based  on  a  -formulation  by  Kinner*  which  has  been  avail¬ 
able  since  1937.  Prandtl  had  pointed  out  that  the  then  new  acceleration  potential 
approach  would  permit  o:.c  to  obtain  exact  solutions  (it  no  longer  required  a  priori 

assumptions  about  the  shape  of  the  trailing  vortex  sheet).  The  circular  wing  in 

2 

incompressible  flow  appeared  to  be  a  suitable  model.  Kinner  -formulated  the 

analysis  and  succeeded  in  determining  overall  lift  and  moment.  However,  the 

Kinner  results  appear  to  diverge  toward  the  wing  tips,  and  the  formulation  appeared 

to  be  too  cximbcrsome  to  be  of  much  practical  value. 

Actually,  solutions  of  the  Kinner  formulation  can  readily  be  calculated  to 

considerable  accuracy  by  solving  a  relatively  simple  linear  system.  The  truncated 

infinite  set  of  numbers  which  one  obtains  converges  to  zero  only  slowly.  The  bulk 

of  this  set  describes  the  solution  for  the  planar  wing  (this  solution  we  have  chosen 

a s  our  reference  solution;  its  chief  characteristic  is  the  wing  tip  singularity). 

With  any  other  solution,  the  part  which  is  of  primary  engineering  interest  is  the 

pressure  distribution,  over  tho  inner  part  of  the  wing;  this  distribution  is  contained 

in  the  void  set,  that  part  of  the  numerical  set  which  is  void  of  the  tip  singularity. 

This  void  set  converges  much  more  rapidly. 

The  analysis  by  which  the  details  of  the  structure  of  the  exact  solution  were 

2-5  2  3 

derived  has  been  described  elsewhere  .  The  leading  terms  were  used  to  cal¬ 
culate  the  initial  shape  of  the  vortex  trail  (and  thus  to  indicate  necessary  corrections 
to  earlier  analyses).  The  complete  solution  is  described  in  Kef.  4;  certain  results 
of  more  general  mathematical  interest  are  presented  in  Kef.  5. 

In  Section  2  of  this  paper,  the  here  relevant  results  of  the  analysis  are 
reported.  In  Section  3,  tho  numerical  methods  are  descr"  1  by  which  wc  have 
calculated  the  reference  solution  to  a  very  high  accuracy;  though  this  accuracy  is 
not  normally  required,  it  servos  to  demonstrate  the  details  of  its  struc  >rc. 
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Further  accurate  solutions,  obtained  with  relative  ease  once  the  reference  solution 
is  known,  are  presented  and  discussed  in  Section  4. 

An  alternative  simplified  approach,  useful  if  one  is  interested  in  the  solution 
ov.’r  the  inner  part  of  the  wing  rather  than  in  the  tip  singularity,  is  presented  in 
Section  5.  Sample  pressure  distributions  are  discussed  in  Section  6. 

The  methods  here  presented  allow  one  to  investigate  a  wide  variety  of 
problems  of  technical  interest.  The  present  paper  is  confined  to  cases  where  the 
given  tlownwash  is  symmetric  and  docs  not  ct.vi.ain  camber.  However,  this  special¬ 
ization  is  not  essential. 

2.  Analytical  Results 

(a)  Pressure  Distribution 

We  describe  the  pressure  difference  p  between  lower  and  upper  wing  surfaces 
by  a  non-dimensional  pressure  function  p.  From  p  the  singularities  in  p  of  orders 
t'j-  at  the  leading  edge  (1.  e. )  and  the  trailing  edgo  (t.  e. )  have  been  removed: 

ps(l-J2)2p/q  (2.1) 

Here  q  is  the  dynamic  pressure  and  the  local  chordwise  coordinate  such  that 
Jl.e.  =  -1’  ?t.e.  =  +1* 

The  given  boundary  condition  is  the  downwash  wV  on  the  wing.  V  is  the 
forward  speed,  w=tand,tA  the  local  incidence.  Consider  the  two  limit  cases  for 
which  solutions  are  known,  the  infinite  uniform  wing  (2-D;  the  two-dimensional  case) 
and  the  slendor  wing  (s.w. ),  specifically  here  the  elliptic  s.w.  in  the  limit  k^O, 
where  k  the  axis  ratio  of  the  ellipse.  The  nature  of  the  chordwise  distribution  p(J-) 
differs  distinctly  between  these  two  limits.  Assume  w  constant,  and  normalize 
to  w=l.  One  then  has  tbe  following  comparison: 


r . . .  . — 

?  4 
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P(J) 

CA 

U 

2-D 

40-J) 

2ir 

-0.5 

s.w. 

* 

r?}  “i 

r 

t  r-ol 

2k 

.tt/4  (2-2) 

where  is  the  local  lift  coefficient,  ^  the  local  center  of  pressure. 

In  2-D,  p(p  is  a  polynomial  of  order  (n+1)  in  j-  if  w(^)  is  a  polynomial  of 
order  n.  Such  smooth  distributions  p(£)  wc  call  2-D  typo  distributions.  In  contrast,, 
the  s.w.  distribution  has  a  kink  at  the  wing  axis  J=0. 

On  a  wing  with  finite  aspect  ratio  and  with  rounded  (parabolic)  tips,  2-D  type 
distributions  nrmst  prevail  over  the  inner  part  of  the  wing  span  but  the  kinked  s.w. 
type  must  exist  at  the  tip  itself.  Indeed,  this  is  confirmed  by  the  exact  solution. 

The  latter  is  derived  from  the  K inner  *  formulation  for  the  distribution  of  the  pressure 
function  p  over  the  circular  wing: 

‘  21  C2Kr2Kco,2ltf  +  ll  C2ulr2A+1,i„(2Ul)^  (2.  3) 

Thus  p  is  built  up  from  elementary  solutions  (none  of  which  is  physically  meaning¬ 
ful  individually)  with  unknown  amplitudes  C^.  Coordinates  and  notations  in  (2.  3) 
are  defined  in  Fig.  1.  *  In  addition  to  (2.  3),  we  have  for  p  the  Kutta  condition 

*  Our  notation  differs  somewhat  from  that  of  Kinner*.  Contained  in  (2.  3)  is  the 
assumption  that  the  down  wash  which  is  given  on  the  wing  is  both  symmetrical  and 
without  camber 

w(|,y)  a  w(y)  •••  w(-y) 

Kinner*  shows  that,  split  off  in  a  suitable  manner,  the  camber  part  p  of  a  general 

c 

pressure  function  is  trivial,  and  zero  all  along  the  wing  edge.  Also,  the  analysis 
can  readily  be  rewritten  for  the  case  that  w(y)  is  antisymmetric. 


»  1 

1 

'  1 
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Pt.e.  5  P*1’  -T</5<0)  5  0 


(2.3a) 


Integrating  (2.  3)  chordwise,  we  obtain  the  local  lift  coefficient 


oo 


Cx(y)  =  (2T/y)  ^  c2i!P2K(y) 

K=0 


(2.4a) 


and  the  local,  moment  coefficient 

bW^(y)  =  •?!  (2»Hte2ulP2A( y) 


(2.4b) 


The  P 2^  are  the  Legendre  polynomials,  the  their  derivatives  (reap,  the  ultra- 
spherical  or  Gcgenbauer  polynomials) 

«  _  1  d 


'  21  "  (X  +  T)  (2^+ 1 )  dy  P2Xd  l(y) 

Both  are  symmetrical  and  orthogonal,  normalized  such  that 

p2jt(il)=  Vtl}  =  1 

Integrating  also  spanwise,  one  obtains  the  vector  of  the  total  lift 

CL  •  «C0  !  L  ”  -C/3C0 


(2.4c) 


(2.  4d) 


( 2.  4e) 


It  is  given  by  the  first  two  amplitudes  alone. 

Before  turning  to  the  analytical  results  which  arise  from  combining  (2.  3) 
and  (2.3a),  let  us  review  two  graphical  presentations  of  representative  results  and 
relate  them  to  (2.2).  These  results,  Figs.  2  and  3,  are  for  our  reference  case, 
the  planar  wing  normalized  to  w(y)  -  1. 

Shown  in  Fig.  2  is  a  relief  of  p  over  the  right  half  wing.  Noting  the  result 
both  of  lifting  line  theory  (large  span)  and  of  s.  w.  theory  (small  span)  that  the 
lift  distribution  over  planar  elliptic  wings  is  an  ellipse,  and  assuming  that  this 
would  hold  also  for  the  circular  wing,  wt  would  La-re  =  const.  =  Cj  .  In  Fig.  2, 

CT  si.  790..  .  Using  rounded-off  numbers,  we  would  thus  expect  p.  =4Ct/TT=2.28 
from  2-D  in  (2).  The  actual  p.  curve  in  Fig.  2  lies  somewhat  higher;  this  is 
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due,  mainly,  to  the  fact  that  the  chordwisa  distributions  p(|)  are  not  linear 
(as  in  (2.2))  but  are  concave  curves.  However,  they  are  2-D  type  curves  as 
expected,  over  the  inner  part  of  the  wing  span. 

For  the  wing  tip  we  should  have  Pj  e  =  Pj,  =  20^  -  3.  58  from  s.  w.  in 
(2.2).  The  actual  value  in  Fig.  2  is  =  3. 19.  That  this  value  is  distinctly 
smaller  than  predicted  is  not  a  failure  of  s.  w. .  Rather,  it  turns  out  that  (as  was 
already  likely  from  Ref.  1 )  the  correct  lift  distribution  over  planar  elliptic  wings 
is  not  an  ellipse. 

Additional  information  about  the  transition  to  the  tip  is  provided  by  the 
curve  for  fa,  Fig.  3.  Over  most  of  the  span,  is  not  much  forward  of  the  {--chord 
point  £  =  -0.  5.  The  transition  to  the  tip  value  -T/  4,  see  (2.  2),  occurs  rapidly 
and  with  a  vertical  tangent. 


The  singularity  at  the  tip  which  is  illustrated  by  Fig.  3  is  not  apparent  in 
the  curve  p.  of  Fig.  2,  but  a  singularity  in  p  occurs  as  the  tip  is  approached 
along  the  wing  axis  £=0. 

A  measure  for  the  amplitude  of  the  tip  singularity  in  p  is  the  tip  value  p^,. 
For  it  one  finds,  considering  that  due  to  (2.  3a)  the  two  sums  in  (2.  3)  must  equal 
each  other  at  the  1 .  e. , 

co 

p„,  5  lim  P(l,^)  =  8^,(2A+1)C?.  (2.5 

1  /}•»+  0  A  =  0  ' 

One  also  shows  that,  moving  toward  the  tip  along  any  half-ellipse  £=const.j  one 
reaches 


lim  p 
y-r  1 


(|=const. ) 


(2.5a) 


This  is  in  agreement  with  the  s.w.  distribution  p(£)  in  (2.  2).  According  to  (2.  7), 
the  p  curves  of  Fig.  2  over  the  half-ellipses  ^  ~  -0.  8,  -0.  6,  , .  reach  the  tip 
vertical  at  equidistant  points  p  =  0.  8p^lf  0.  6p,j,,  . .  .  All  those  for  ^0  reach  the 
same  end  point  p=0. 
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Due  to  (2. 4b),  the  malr  be  called  the  moment  amplitudes;  according 

to  (2. 5),  they  determine  p^,.  For  the  lift  amplitudes  we  have  the  aero  sum 
condition 


(2.6) 


Utilizing  the  orthogonality  of  the  circular  functions,  one  finds  from  the 
equality  along  the  1.  e.  of  the  two  sums  in  (2. 3)  that  the  two  sets  of  amplitudes  are 
fully  interdependent  and  equivalent.  Introducing  the  abbreviation  I=(A+£),  we  have 

m  A C„ 


Tcatl  -  2  X 


*=0A2-*2 


(all  A) 


(2.  7i) 


ire,..  =  2  t  -AP- 


(2.7b) 


The  transformation  (2. 7o)  leads  to  a  set  which  automatically  fulfills  (2. 6). 
As  an  example,  set  Cj=l;  for  A/0.  Then 

'i?C0  =  2  ;%,0=V-.A»-’  (2-8) 

One  readily  shows  that  this  set  fulfills  (2.6). 

The  set  (2.  8)  is  one  of  a  sequence  of  possible  elementary  sets  which  one 
might  use  Ic  build  up  the  complete  solution  for  a  given  dowmvash  w(y).  Choose  an 
integer  A*  and  let 


C2A*+1  ^  °  !  C2Hl=0for^A* 


(2.9) 


The  pressure  distribution  which  results  is  characterized  by  its  1.  e.  distribution 


and  is  thus  a  spanwise  wavy  distribution.  We  will  further  discuss  the  elements 
A*  in  Section  6. 
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The  elements  A*  are  not  an  optimum  choice  because  one  reqxiires  an 
infinite  number  of  them  to  produce  a  downwash  w(y)  which  is  regular  at  the  wing 
tip.  We  generally  denote  by  En  any  elementary  set  which  we  use  to  build  up  the 
solving  set  C,,^,  and  we  always  require  that  E^  individually  fulfills  the  zero  sum 
condition  (2.6).  From  (2.8)  arises  the  suggestion  to  consider  the  sets 

Eo*$(r)  ;  -  K  r  W/ 0)  (2.  10) 

Here  J denotes  the  Riemann  function,  see  c.g.  Table  23.3  of  Ref.  6.  We  call  the 

sets  Ek  the  rational  sets;  we  will  also  have  to  xxse  logarithmic  sols. 

2 

According  to  (2.7b),  the  set  E^  is  the  leading  set.  Thus 

a. 


'2»< 


? 

-  +  Ry 

1C2  K 


(k/0) 


(2.11) 


In  (2.  11),  *s  rjl  unknown  amplitude,  and  the  remainder  R^  is  defined  by  the 
-2  • 

condition  R  «o(R  ).  One  can  show  that  R^  does  not  contribute  to  the  tip  value 
p^,  and  thus  does  not  contribute  to  first  order  to  the  i*p  singularity.  Hence  we  are 
mainly  interested  in  the  elementary  pressure  function  p  which  arises  from  inserting 

M 

2  2  2 
Ek  into  (?..  7a)  to  determine  its  corresponding  set  (  j,n«d  then  inserting  both 

2 

and  +  j  into  (2.3).  For  we  have  analytical  results  for  both  the  1.  e.  and  the 

wing  axis  fi  -Os 

Pz.i.c.  *  <7-  f ' 

-  4  -  r<y2>3/y  (212> 

"  4y  £  21og(l/ y)  +  1  +  y2log(l/  y)  +  . .  .1 


The  function  f  here  is  the  dilogaritlun,  see,  e.g. ,  27.  7  of  Ref.  6.  The  two  functions 
(2.12)  exhibit  the  characteristic  properties  of  p  in  Fig.  2:  a  smooth  curve  along 
the  1.  e. ,  a  vortical  tangent  at  the  tip  in  P.txjs- 
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From  (2. 11)  and  (2.  )2)  follows  Pj,  =  Sd'a*,.  From  (2.2)  and  (2.  5a)  follows 
^  =  20^(1).  Hence 

Ce(l)  =  Px/2  =  4,,a2  (2.13) 

(b)  Structure  of  the  Solution 

We  describe  the  given  dowmvash  \v(y)  by  a  set  of  downwash  coefficients  w 


co 


w(y )  =  21  (s+l)(2sil)wfiP2s(y) 


6  =  0 


(2.  14) 


where  the  {^(y)  arc  Gcgeubauer  polynomials  °f  (2.4c)  and  are  subject  to  the 
normalization  (2.4d).  Reversing  (2.  14)  we  have 


\v 


= _(4«+3)  J  {  ]  ^2b+l(y)fljr 


(2. 14a) 


o  o 

For  example,  in  our  reference  problem 


wo =  1  ;  we/o  =  0 

The  relation  between  the  sets  w  and  is* 

s  ZK 


where 


w  = 

s 


S,(C,  >  =  i  'k.  P-T  -  l«-CZ)t  <2. 15a) 

1  ^HallsfH+s)  (s+i)(Ms+J)J  Zlt 


=  *C2»+l  +  S)<C2>i> 


(2. 14b) 


(2.  15) 


The  set  C2s+  j  is  given  by  (2,  7a),  with  A  replaced  by  s.  The  operator  S^  is  of 
higher  order;  the  abbreviation  i=s+|  corresponds  to  l. 

One  obtains  information  about  the  structure  of  the  set  by  observing 
that  in  all  technically  meaningful  cases  the  downwash  w  is  regular  at  the  tip.  From 
this  follows  that  the  set  w  must  converge  more  rapidly  than  any  power  of  (l/s). 


*  This  form  is  considerably  more  convenient  analytically  than  the  original  pre¬ 
sentation,  Eq.  (CO)  of  Ref.  1.  . 
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To  this  requirement  we  refer  as  the  tip  condition  for  the  set  wg. 
it  *  r 

The  C2s.|  j  which  arise  from  the  1'^  and  which  thereby  occur  in  (2.  15)  are 

.3 

generally  0(s  );  hence  there  have  to  be  certain  relations  between  their  amplitudes 

such  that  terms  of  this  order  (as  well  as  of  higher  orders)  cancel  out.  The 
investigation  becomes  fairly  intricate  because  of  the  occurrence  of  logarithmic 
sets  in  addition  to  rational  sets;  on  the  other  hand,  it  is  those  logarithmic  sets 
which  enforce  certain  useful  general  relations.  The  result4  we  write  here  in  the  form 


-C 


2* 


lK(4>(+ 


31°gV 
K(4>M-1)  16S2K4J 


a2  + 


(2. 16) 


Abbreviating,  we  write  this  as 

"C2)t  =  +  RK  (2.16a) 


and  describe  it  by  saying  that  is  composed  of  a  general  set  GK,  common  to  all 
solutions,  and  a  remainder  set  R^.  For  the  unknown  amplitude  C^(J.)  of  we  have 
(2.6)  with  (2.  13). 

Between  (2.  11)  and  (2.  16),  certain  elementary  sets  have  been  transferred, 

by  means  of  the  tip  condition  for  the  set  ws,  from  the  remainder  R^  of  (2. 1 1)  into 

the  general  set  G^.  On  principle,  one  could  continue  this  process.  However,  this 

would  serve  no  useful  purpose.  Namely,  the  remainder  of  (2. 16)  (which  differs 

of  course  from  that  of  (2. 11))  contains,  in  addition  to  further  general  sots,  also 

specific  sets  which  describe  the  solution  for  the  specific  given  downwash  w(y).  The 

-4  -4 

orders  of  magnitude  of  the  two  leading  specific  sets  are  0(K.  logx)  and  0(K.  ). 

There  is  no  practical  advantage  in  knowing  general  sets  which  are  of  higher  order 

-5  2 

than  these  specific  sets  (the  next  general  set  is  presumably  0(K.  log  w)). 

We  will  require  also  some  information  about  the  structure  of  the  sets  Cg^j. 
If  we  rewrite  (2.  16)  in  the  alternative  form 


•3)og2»t 

i6?rV . 


a2  +  Rrt. 


(2. 16b) 


WKW9Ntf0WBSi 
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we  have  again  R  =  0(K  JogK.)  anti  have  for  large  l 


it  c. 


(2. 16c) 


:2.Wi  =  {[l  +  I-A  -  log  2]  «2-2|,*R1J/a-» 

*3a2log/i/l6?i4  +  °(1/A4) 

The  set  is  a  logarithmic  set  which  is  related  to  tliu  V' (digamma)  function;  it  is 


defined  by 


V  1  +  T  +  5  +  •  •  •  +  2j[TT  *  i  Ifa+i » -  V<1)] 


-  |logA 


(2. 16d) 


(c)  Span  Loading  and  Tip  Distribution 

In  order  to  translate  (2.  16)  resp.  (2. 16b)  into  a  formula  for  the  lift  coeffi¬ 
cient  Cj(y),  we  need,  according  to  (2.4a),  the  functions 


e(y)  - 


EMP2K<V' 


(2.17) 


which  are  created  by  the  various  elementary  sets  in  (2.  16).  We  have  analytical 

5 

results  only  for  the  two  leading  sets 


e2<y> s  2y  +  >y?‘iog(2/y)  +  c%2) 

e3(y)=  y2log(2/y)  +  0(y?) 


(2. 17a) 


Using  these  in  (2.4a  with  (2. 16),  we  obtain 

C*(y)  -  [l  +  ;i/8)ylog(2/y)]  0^(1)  +  0(y)  (2- 18> 

A  useful  application  of  this  formula  is  illustrated  in  Fig.  4..  Numerical 
summation  of  (2.4a)  becomes  cumbersome  near  the  tip  because  of  the  denominator  y. 
Shown  in  Fig.  4  are,  near  the  tip  and  to  a  large  sca?e,  the  approximate  curves  which 
have  been  calculated  truncating  (2.  4a)  after  N  =  20,  25. .  40  term;'..  The  last  two 
curves  have  converged  well  enough  io  the  left,  that  is,  inside  y=t0.  96,  but  all  curves 
show  increasing  waviness  toward  the  right  and  turn  up  toward  l-ao  at  the  tip.  This 


WWJBWMiW 
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remains  true  for  all  finite  N.  Also  shown  is  the  asymptotic  curve,  the  term  with 
C^(l)  of  (2, 18).  To  obtain  a  fairly  accurate  curve  (^(y),  one  has  to  match  the 
asymptotic  curve  with  the  curve  N=40  near  y=0.  98  and  has  to  smooth  out  the  latter 
curve  between  y=0.  96  and  0.  98. 

We  insert  here  an  illustration,  Fig.  5,  of  a  similar  matching  process  for 
the  cliordwise  distributions  p(p  near  the  wing  tip.  The  curves  shown  represent  the 
result  of  using  in  (2.3)  the  first  N-40  amplitudes  of  each  set,  and  +  j  of  our 
reference  solution.  Again  convergence  deteriorates  as  y-M  because  of  the 
denominator  y.  The  largest  error  due  to  truncation  occurs  along  the  wing  axis. 

Here  the  points  due  to  p^  alone,  ace  (2. 12),  are  shown  as  circles. 

When  y=0.  980,  curve  and  point  almost  coincide.  The  truncation  error  is 
still  small,  and  the  leading  contribution  8 till  almost  equals  the  total.  As  y  is 

increased,  the  latter  becomes  increasingly  more  correct  while  the  truncation  errors 
increase.  The  slight  wavinoss  in  the  curve  for  y=0.  995  is  presumably  not  genuine 
but  due  to  truncation  errors. 

Both  Fig.  3  and  Fig.  5  show  that  the  transition  from  the  kinked  tip  distri¬ 
bution  to  the  smooth  2-D  type  distributions  occurs  rather  rapidly. 

With  Fig.  6  we  return  to  span  loadings.  The  reference  planar  solution  is 

shown,  and  all  other  curves  are  normalized  to  its  tip  value  (1).  Illustrated  is 

2 

the  role  of  the  leading  terms  in  (2. 16).  The  curve  (y)  represents  the  contribution 
o 

of  E‘;  in  it,  the  factor  (l/S)  of  (2. 18)  is  replaced  by  (1/ 4),  see  (2. 17a),  and 

accordingly  it  cuts  over  the  correct  curve  near  the  tip. 

2 

The  downwash  which  belongs  to  is  shown  in  Fig.  7,  denoted  as  It  is 
negative  infinite  at  the  tip;  that  of  the  elliptic  lift  distribution,  not  shown  in  Fig.  7, 
is  positive  infinite  at  the  tip. 

In  Figs.  6  and  7,  basic  refers  to  the  simplest  set  which  produces  a 
finite  value  w(l): 


n  ^4 1 

*C2K,  basic  ~(.K2  "  4^3  +  k4  i 


a. 


(2.19) 


k 
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The  first  two  terms  on  the  right  cancel  between  them  a  term  0(s  log  s)  in  w  (and 

s 

produce  the  correct  factor  (1/ 8)  of  (2.  18)).  Their  combined  downwash  is  shown 

_3 

in  Fig.  7.  The  third  term  in  (2.19)  is  added  to  cancel  also  the  term  0(s  )  in  w  . 

6 

* 

The  amplitude  a.  which  achieves  this  is 


1®+U?L. 

a4  r  fSJUH-  4f(3) 


0.3741083... 


(2.19a) 


The  basic  downwash  curve  in  Fig.  7  is  seen  to  be  almof:t  a  horizontal  line,  but  it 
has  an  infinite  tangent  at  the  tip.  The  curve  basic  in  Fig.  6  goes  asymptotically 
into  the  planar  curve  at  the  tip. 

The  point  marked  collocation  analysis  in  Fig.  6  is  determined  from  numcri- 
7 

cal  results  ,  calculated  by  means  of  an  unusually  careful  spanwise  integration. 

The  collocation  results  coincide  with  the  correct  planar  curve  over  the  inner  part 
of  the  wing  span.  It  deviates  at  the  tip  by  5.4%;  this  deviation  arises  because 
collocation  analyses  assume  an  elliptic  type  lift  distribution  and  thus  a  finite  slope 
in  the  C^(y)  curve  at  a  parabolic  wing  tip. 

The  pressure  formula  (2.  3)  and  the  result  (2. 16)  are  specific  to  the  circular 
wing,  but  insofar  as  (2. 16)  can  be  translated  into  the  analytical  form  (2.  18),  it 
can  be  generalized  to  apply  to  arbitrary  wings  with  parabolic  wing  tips^.  For  such 
wings  the  local  lift  becomes 

X  (y)  -  |[l  +  (^-)*log-p-J  c£(l )  +  R(y)J  c(y)n  (2. 20) 

Here  b  is  the  wing  span,  r  the  wing  tip  radius,  c(y)  the  wing  chord.  The  remainder 
function  R(y)  may  be  expected  to  have  a  finite  tangent  at  the  wing  tip. 


3.  Numerical  Procedure 

For  a  given  downwash  w(y)  one  finds  the  set  w  by  means  of  (2.  14a)  and 

8 

has  in  (2.  15)  an  infinite  linear  system  for  calculating  the  set  C^,  Solving  this 

N 

system  after  truncating  it  to  N  equations,  one  obtains  N  values  which  represents 
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a  truncated  approximation  to  the  correct  set 

Convergence  toward  the  correct  set  slows  down  as  N  is  increased 

because  the  set  itself,  being  a  reciprocal  progression,  converges  with  in- 
6H 

creasing  slowness  as  K->oo.  To  improve  the  convergence,  one  calculates  first 
the  set  C^+j.  For  this  purpose  (2. 15)  is  rewritten  as 

w„ "  IT  _  L  I  ~  '  TT  2-  +1  (5*1) 


'V  Cxc2x+isirc28+i  c^Cz>+1  (3,  l) 

**  s 

The  matrix  c®  is  almost  the  unit  matrix;  the  numbers  c ^  (wh-'ch  correspond  to  the 
operator  Sj  in  (2. 15))  are  small.  A  short  numerical  table  follows: 


0 

1 

2 

3 

1 

0 

0.9112 

-0.0470 

-0.0315 

-0.0236 

i 

-0.0117 

0.9902 

-0.0078 

-0.0064 

2 

-0.0039 

-0.0041 

0.9964 

-0.0032 

3 

-0.0018 

-0.0021 

• 

-0.0021 

0.9981 

The  c®  are  given  by 


Here 


c  =  MS  -  (M-N)L 

® » A 

=  T  -  [4i/(4i+l)JS 


if  A/s 
if  A=s 


M  =  i/(X2-(s+l)2)  ;  N  =  i/(JZ.|2) 

MM  )  =  21^  +  (1/ 2X)  -  (l/2s)  ►  2Lg 

S  *  21°e2  ‘  ^  (27-T)  '  2(2s+  l){srrj 


IMS2/4)-2-l_.-> 

vxi  (2*>-i r  (2s+ir 


(3.2) 


(3.2a) 


The  set  L  is  defined  by  (2.  l6d). 
s 


15 


Having  calculated  the  set  C?A+J  using  a  truncated  system  (3. 1),  one  uses 

N 

(2.7b)  to  calculate  the  corresponding  truncated  set  C.,^.  * 

The  complete  procedure  is  easily  programmed  for  digital  computation**. 

No  large  number  N  is  required  to  obtain  engineering  accuracy.  Ono  obtains  an 
accuracy  of  about  10  with  N-40  in  the  reference  problem  (2. 14b);  of  course,  the 
accuracy  goes  down,  or  N  has  to  be  increased,  if  higher  coefficients  w  occur.  With 
Ns40,  one  has  more  than  sufficient  accuracy  for  moat  purposes  over  the  inner  part 
of  the  wing.  If  one  wants  to  match  near  the  tip,  as  in  Figs.  4,  5,  one  needs  also 
the  amplitude  the  tip  singularity.  If  no  unusual  accuracy  is  required,  it  is 
sufficient  to  assume  that  R^=0  in  the  highest  C  ^  which  have  been  calculated;  then 
is  obtained  directly  from  (2. 16). 

4.  Calculation  of  the  Reference  Solution 
(a)  Approach 

Wc  have  determined  the  reference  solution  plus  three  further  solutions  to 
very  high  accuracy  in  order  to  definitely  determine  and  demonstrate  the  details  of 
the  structure  of  the  solutions  (and  also  the  powerfulness  of  the  available  methods). 
We  used  up  to  N=60  equations,  used  12  decimals  in  the  results  to  extrapolate  to 
N=oo,  lost  2  decimals  in  the  extrapolations  and  performed  the  further  calculations 
with  10  decimals.  The  resulting  sets  Ate  listed  in  Tables  I  and  II,  rounded  to 
8  decimals. 

The  C^K  which  arc  calculated  directly  from  (2. 15)  decrease  as  N  is  increased; 
those  calculated  via  (3. 1)  and  (2.  7b)  increase  with  N.  The  two  different  results 
thus  establish  reliable  (though  not  overly  narrow)  bounds. 

**  A  program  for  the  UNIVAC  1108  which  delivers  both  C^+j  and  C^K,  with  N 
a  variable,  can  be  obtained  from  the  author. 


It  is  easy  to  extrapolate  to  N*co  results  which  converge  like  geometric 
series.  In  the  present  problem  behavior  like  reciprocal  powers  of  N  was  more 
likely  (sec  Appendix  £  of  Ref.  8)  and  was  in  fact  indicated  by  the  numerical  results. 
For  the  +  y  we  «8®d  successfully  the  following  principle:  Assume 

f(N)  =  f(en)  *  A/Nr  (4.1) 


where  A  and  r  are  unknown  constants.  Then  to  second  order 

2A(N)  =  Jf(N-tl)  -  f(N-l)  =  2rA/Nr+1 


Thus 


from  which  one  finds  r,  and 

f(oo)  -  f(N)  +  NA(N)/  r 

The  point  now  is  that,  while  we  do  not  know  if  (4. 1)  correctly  describes  the  dependency 
N 

of  on  N,  the  same  procedure  works  fairly  well  if  the  term  on  the  right  should 

bt  AlogN/Nr,  for  example.  In  fact,  this  seemed  to  be  the  case,  with  r  of  the  order 
r«3.  By  performing  the  evaluation  wilh  several  values  M,  one  obtains  information 
about  the  validity  of  the  assumption  made  and  about  the  accuracy  of  the  resulting 
value  f(oe>). 

The  set  C ^  required  a  modified  approach  to  be  described  below.  The  final 
set  was  used  to  determine  the  ampoilude  C^(l)  of  the  tip  singularity. 

(b)  Outer  Solution 

We  call  the  amplitude  Cg(l)  the  outer  solution;  it  has  to  be  determined  by  a 
matching  process.  The  difficulty  of  matching  in  term3  of  C^(y)  is  illustrated  by 
Fig.  4.  A  more  definite  approach  is  provided  by  the  sum  condition  (2. 5),  which  we 
write  here  in  the  form 


C^l)  n  S(N)  i  a(N)  s  4  |(21+1)C2mi  +  4  £  (2*+ 1)0^  (4.2) 


One  visualises  that  one  could  calculate  S(N)  for  a  number  of  values  N,  plot  these 
over  u=(l/N),  and  extrapolate  to  the  origin  u=0.  The  end  value  of  the  extrapolation 
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curve  would  then  be  C^(l).  However,  we  know  from  (2.  16c)  that  the  leading 
elements  in  the  set  to  be  sximmed  have  the  form 


<2X+1)C2X+ 


1 


DQlogX  +  Dj  D2logX+  D3 


<2Xh  1)' 


(2X+1T 


+  •  •  • 


(4.3) 


Therefore  A(N)  behaves  like  ulogu;  S(N),  plotted  over  u,  will  have  a  vertical 
tangent,  and  extrapolation  is  unreliable. 

To  eliminate  this  difficulty,  replace  (4.2)  by 


with 


C/(1)  =  S*(N)+  A*(N) 

S*(N)  »  S(N)  +  4N(2N+1)C2N+1 


Then,  using  formula  (A4)  of  the  Appendix 


(4.4) 


(4. 4a) 


A  (N)  = 


4N+3 


(2N+1) 


r  D0  -  [(log  N  -  i)D2  1  D3] 


_ 2N 

(2N+1)3 


(4.4b) 


The  coefficient  Dj  has  disappeared  from  the  leading  terms,  and  the  curve  S#(N) 

has  the  form  C^(l)  -  0(u).  This  curve  was  used  in  Ref.  2.  It  is  shown  at  the  far 

-5 

left  of  Fig.  8;  the  error  in  its  end  value  C^(l)  was  about  10  . 

The  marked  improvement  which  is  illustrated  by  the  two  main  curves  of 

4 

Fig.  8  is  due  to  matching  with  new  analytical  results  .  According  to  (2. 16c) 

D0  =  C^D/T2  J  D2  =  -3C>(1)/41?  (4.5) 


Inserting  at  first  only  Dq  into  (4. 4b),  we  replace  (4. 4)  by 

'(N)/[l  -  -y  4'N+3 
l  1r  (2N+ 1 ) 


<*<»  *  S' 


4m3  V2j.OBN.-b, 

«  J  i'iVT  i  1 


(2N+1 ) 


(4.6) 


with  constants  a  and  b.  The  first  term  on  the  right  we  denote  by  S*/[],  Its  curve 

in  Fig.  8  has  to  have  a  horizontal  tangent  at  its  origin  C^(i). 

The  lower  curve  in  Fig.  8  has  been  constructed  using  the  original  sets 
N 

^2A  +  1  ra'k°r  than  the  extrapolated  set;  this  curve  must  reach  the  correct  origin 
C^(l)  but  we  do  not  know  if  it  has  to  have  a  horizontal  tangent. 
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The  upper  curve  was  calculated  using  the  final  set  C2A+1»  uPt°N=40. 

2 

(Accuracy  deteriorates  when  N  becomes  large  because  of  the  factor  N  in  (4.4a). ) 
Disregarding  all  higher  order  terms  in  (4.6),  the  three  constants  0^(1),  a  and  b 
were  determined  such  that  the  resulting  curve  passed  through  all  points  (l/N).  The 
agreement  between  this  curve  and  the  calculated  points  extended  down  to  N=5,  far 
to  the  right  of  Fig.  8.  The  origin  of  this  curve  is  C^l)  =  1. 5930904. .  . 

This  agreement,  while  it  tends  to  be  reassuring,  must  bo  somewhat  coinci¬ 
dental,  however;  it  is  highly  unlikely  that  the  higher  terms  in  {4. 6)  should  be 
negligible  down  to  N=5.  Indeed,  we  should  correctly  have  see  (4.  5).  This 

condition  was  not  quite  fulfilled,  and  an  adjustment  had  to  be  made,  allowing  for 
small  but  finite  higher  ordor  terms  in  (4.6).  The  adjustment  led  to  the  new  value 
0^(1)  =  1.  5930884. .  .  Both  new  values  are  marked  in  Fig.  8.  The  smallness  of 
the  difference  between  them  may  be  taken  as  an  indication  of  the  degree  of  accuracy 
which  had  been  achieved. 


In  Table  I  the  value  ag  =  0.  1267740  is  used;  to  it  corresponds  C^(l)  =  1.5930890. 

In  the  context  of  Table  I,  the  number  has  to  be  treated  as  an  exact  number;  however, 

•  • 

tho  listed  value  is  not  claimed  to  be  correct  beyond  the  accuracy  range  which  is 
indicated  by  Fig.  8. 

(c)  Listing  of  the  Solution 


The  reference  set  1®  listed  3n  Table  I  where  it  is  denoted  by  C ^ 

We  discuss  first  how  this  table  was  calculated. 

N- 1  N 

Two  successive  sets  and  C-  ^  differ  because  of  two  different  reasons. 

N-l  N  N 

One,  the  sets  and  Cg^+j  differ.  Two,  *‘.e  new  end  term  enters  in  (2.7b) 

N 

when  i®  calculated.  The  second  cause  has  by  far  the  larger  effect.  Therefore, 
we  wrote  for  N  large 


rN-i  1 


-  K 


(times  jr  if  ^-0) 


(4.7) 


The  number  €  is  here  introduced  to  take  care  of  reason  one. 
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We  found  that  4  changed  little  with  H  and  less  with  N.  If  one  assumes  that, 
as  N  increases,  e  remains  constant  and  C^N+l  behaves  like  logA/  A  (compare  (4,  3); 
the  constant  is  small  in  the  case  of  the  reference  solution)  one  can  sum  the 
differences  (4.  7)  to  infinity,  and  thus  obtains  C^.  Doing  this  for  both  N=40  and 
N=60,  we  found  that  our  results  were  the  same  within  our  accuracy  requirements 
(the  agreement  began  to  deterioriatc  beyond  K=25  since  k approached  N=40,  but 
results  beyond  »t=25  were  not  required). 

Table  I  lists  C ^  up  to  h=15,  and  lists  the  remainder  which  was  cal¬ 
culated  by  means  of  the  listed  value  a The  remainder  is  broken  down  further 
in  the  manner  shown  (J*K  is  again  the  logrithmic  set  defined  by  (2.  I6d)),  and  a 
final  remainder  R*  is  listed  in  the  last  column.  This  column  is  zero  from  K=15 
onward. 

The  constants  and  were  determined  by  means  of  the  three  conditions 
e  Pa  = 

(a)  Rp  =  0  ;  (b)£.R^  =  0  ;  (c)R^=0  for  »f>p  (4.8) 

Here  p  is  a  number  which  has  to  be  determined  by  trial  and  error.  One  chooses 
p  and  determines  and  c^,  using  conditions  (a)  and  (b).  The  zero  sum  condition 
(2.6)  must  apply  to  the  set  R^  (that  is,  if  both  analysis  and  numerical  work  are 
without  error)  because  it  applies  individually  to  all  other  sets  in  C2k*  Hence  the 
sum  of  all  R^  from  p+1  to  co  must  be  zero  when  (a)  is  fulfilled.  However,  this 
does  not  insure  that  (c)  is  fulfilled.  Rather,  if  p  was  chosen  too  small,  the  RK 
beyond  p  will  describe  a  wavy  curve  about  the  zero  line. 

On  the  other  hand,  significant  decimals  are  lost,  and  the  constants  b^  and 
c ^  will  be  ill-defined,  if  p  is  chosen  too  large.  In  fact,  it  is  not  possible  to 
determine  these  constants  very  accurately;  though  seven  decimals  are  listed  in 
Table  I,  the  last  decimals  should  not  be  considered  reliable.  A  small  change  in 
b^  and  will  affect  the  first  few  listed  values  R^;  however,  since  these  first  few 
values  almost  cancel  each  other,  the  changes  will  cancel  out  entirely,  and  the 
validity  of  (4.  8)  will  not  be  affected, 
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On  the  other  hand,  the  fact  that  li  was  possible  to  find  p,  bg  and  such 
that  (4.8  )  was  fulfilled  to  10  decimals  (allowing  for  rounding-off  errors)  and  that  the 
set  R^  converged  smoothly  and  rapidly:  this  fact  can  be  considered  an  overall  con¬ 
firmation  of  analysis  and  numerical  work.  (The  value  p  which  we  used  was  p=20. ) 

The  set  q  converges  slowly;  written  out  to  eight  decimals,  it  would 
extend  to  about  )i  =  5,  000.  Table  I  demonstrates  that,  knowing  the  structure  of  the 
solution,  one  can  store  it  in  form  of  the  short  set  R^plus  the  three  numbers 
b^  and  c..  (These  numbers  have  to  be  treated  as  exact  numbers,  of  course,  when 
one  uses  them  to  reconstruct  the  set  (y  )  The  stored  numbers  contain  all  the 
specific  information  about  the  specific  solution;  the  bulk  of  the  numbers  in  the 

9 

infinite  set  C ^  q  describes  the  general  tip  singularity. 

We  close  thissection  on  the  reference  solution  by  noting  also  the  values 
q  and  ^  q  which  arise  by  means  of  (2.  4e).  They  describe  the  vector  of  total 
lift  for  the  planar  circular  wing.  A  fair  number  of  approximate  values  are  given 
in  the  literature,  and  some  of  the  differences  between  these  are  sizeable.  Our 
results 

CT  1.79002303  ;  -J;T  n  *  0. 52085758 

L,o  n.,0  (4 .9) 

should  be  reliable  to  less  than  one  unit  in  the  last  decimal  given. 

5.  Calculation  of  Additional  Solutions 
(a)  Procedure 

To  engineering  accuracy,  solutions  for  a  given  downwash  w(y)  are  obtained 
directly  from  (3. 1).  If  one  wants  very  high  accuracy  in  order  to  treat  more 
details  of  the  solution,  one  can  make  use  of  the  given  reference  solution  q  to 
simplify  the  numerical  procedure.  Write  instead  of  (2. 16) 

-C2*  =  -XC2,(,  0  +  H*  (5.1) 

Here  C2^  q  takes  the  place  of  the  general  set  G^.  In  order  to  insure  that  the  set 
plays  the  role  of  the  remainder  R^,  we  have  to  determine  the  factor  X  such  that 


Cl* 
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*2  =  Xa2, 0  (5* la) 

The  leading  set  in  HR  is  then  again  0(k"^1o£k). 

We  write  Hh  rather  than  RK  in  (5. 1)  to  indicate  that  there  is  a  difference. 
The  set  RK  is  a  remainder;  it  contains  the  specific  sets  which  describe  the  specific 
solution  C^,  but  it  contains  also  those  general  sets  of  higher  order  which  are  not 
listed  in  the  general  part  of  (2. 16).  On  the  other  hand,  the  set  HK  is  void  by 
definition  of  all  general  sets  which  might  be  transferable  from  into  G^.  The 
void  set  Hk  is  that  specific  signature  of  the  given  downwash  w(y,  which  distinguishes 

C2nIroraCZ«,0- 

The  advantage  of  replacing  (2. 16)  by  (5. 1)  is  that  it  is  much  easier  to 

determine  the  factor  X  than  to  determine  the  amplitude  a^.  Our  numerical  pro- 

N  N 

cedurc  corresponded  essentially  to  (4.7).  We  calculated  sets  and  ^  up 

to  N=60  and  determined  a  preliminary  value  X  by  assuming  in  (5. 1 )  that  HR  was 

N 

negligible  for  tffcN.  This  allowed  one  to  calculate  HK  and  to  use  it  in  the  place  of 

N  N 

C2R  on  the  left  of  (4. 7).  In  the  term  which  then  corresponds  to  ^  on  the  right 

of  (4.  7),  the  logarithmic  term  with  Dq  from  (4.  3)  has  been  eliminated  by  (5.  la), 

and  the  rational  term  with  Dj  (no  longer  small)  has  become  the  leading  term.  This 

simplifies  the  summation,  and  hastens  its  convergence.  The  new  extrapolative 

corrections  were  generally  smaller  than  those  of  the  reference  solution. 

The  results  are  listed  in  Tabic  II  in  a  form  corresponding  to  Table  I: 

_  1 


-C2„.nv  =  V^(b4VC4>+H* 


(5.2) 


Here  n  is  the  number  of  the  solution,  and  the  suffix  v  denotes  its  void  part.  To 
find  the  amplitudes  b^  and  c^  requires  »  process  corresponding  to  (4.8).  A  slight 
complication  arose  from  the  fact  that  thv  calculation  had  started  from  a  preliminary 
value  X.  We  had  to  determine  the  effect  of  a  unit  change  in  X  and  had  to  adjust  X 
by  linear  superposition  such  that  the  final  set  H  converged  properly. 
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(b)  Sample  Results 

If  one  stores  the  solutions  for  a  number  of  linearly  independent  downwash 
distributions  wn(y),  emc  can  compose  the  members  of  a  class  of  solutions  as 
required.  A  logical  choice  for  the  wn(y)  are  the  Gegcnbauer  polynomials: 

wn()r)  =  *2n(y)  (5.3) 

The  corresponding  downwash  coefficients  w  arc  from(2. 14) 

s 

Wn,n  =  [(n+D(2n+l)]"1  ;  w^  =  0  fo:  s/n  (5.3a) 

Our  reference  solution  fits  into  this  system  as  the  case  n=0.  We  calculated  in 
addition  the  solutions  for  n=l,  2  and  3.  These  are  listed  in  Table  II. 

The  downwash  wn(y)  i8  given  in  Table  II  in  the  form 

W».v(y,  =  wn<V) -Xn  (5.4) 

The  complete  solution  C,  is  given  by 

C2tf,  n  ~  XnC2|f,  0  +  C2)f,  nv  (5.5) 

The  void  solution  C-  _ is  obtained  from  Table  II  by  means  of  (5.  2). 

CKi  nv  ' 

In  the  columns  of  Table  II,  there  is  always  a  sign  change  between  K=n  and 
K=n+1.  The  numerical  values  in  the  sets  H*  increase  as  n  increases,  and  conses 
quently  the  listed  sets  become  longer.  (We  used  p=25  in  (3. 10)  in  all  three  cases 
of  Table  II. ) 

The  void  span  loadings  ny(y)  which  belong  to  Table  II  are  shown  in  Fig.  9. 
They  appear  to  have  finite  tangents  at  the  tip;  however,  this  has  not  been  investi¬ 
gated  closely.  An  interesting  observation  of  physical  interest  is  the  following:  one 
can  form  the  ratio  between  two  integrals  over  the  wing,  the  total  lift  and  the  integral 
over  he  downwash  momentum,  This  ratio  turned  out  to  be  almost  invariant  between 
the  reference  solution  and  the  three  individual  void  solutions  of  Fig,  9;  the  maximum 
variation  was  less  than  1%.  It  appears  that,  as  far  as  the  total  lift  is  concerned, 
the  details  of  how  the  downwash  is  distributed  over  the  wing  span  arc  of  negligible 
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importance;  they  do  not  seem  to  affect  the  outer  flow  which  determines  the  total  lift. 

The  ratio  varied  somewhat  more  when  the  complete  solutions,  Fig.  10, 
were  compared.  Small  differences  between  large  numbers  occurred  in  the  super¬ 
position,  and  the  lift-momentum  ratio  was  found  to  vary  from  0.  895  for  n=0  to 
0. 821  for  n=3. 

Shown  in  Fig.  10  are  both  the  downwash  distributions  wfi(y)  and  their  span 
loadings  C.  (y).  All  the  latter  curves  form  a  downward  hook  on  approaching  the 
wing  tip,  reaching  a  vertical  tangent  at  the  tip  itself  .  However,  on  the  scale  of 
Fig.  10  this  hook  is  no  longer  visible  when  n=2;  even  though  the  bulk  of  the  set 
still  represents  the  tip  singularity,  the  latter  now  plays  a  minor  role  in  the  span 
loading  curve. 

The  scale  ratio  between  w  and  Ca  in  Fig.  10  is  so  chosen  that  v/rt  and 

n  *,n  U 

o  have  roughly  the  same  overall  magnitude.  For  n /  0,  the  curves  are  wavy, 
and  the  wave  amplitudes  in  C .  are  seen  to  be  distinctly  smaller  than  the  cor- 

^  i  XI 

responding  amplitudes  in  wn<  Each  amplitude  ratio  corresponds  roughly  to  the 
ratio  between  half -wavelength  and  wing  span.  This  is  what  one  would  expect  from 
either  lifting  line  theory  or  slender  wing  theory;  according  to  either,  the  lift 
coefficient  becomes  proportional  to  the  aspect  ratio  when  the  latter  becomes  small 
vCnough,  in  agreement  with  our  observation,  supposing  only  that  one  can  interpret  the 
length  of  a  half-wave  as  an  effective  wing  span.  That  one  can  do  this  (roughly)  is 
not  entirely  unexpected,  but  it  is  of  interest  to  see  the  point  verified. 

Lifting  line  theory  overpredicts  (at  small  aspect  ratios)  by  the  factor  2. 

This  is  illustrated  for  n=3  in  Fig.  11.  Thus  Fig.  11  implies  that  the  correct  result 
corresponds  closely  to  what  we  would  expect  from  slender  wing  theory  for  n=3. 

On  the  other  hand,  there  is  no  slender  wing  type  analysis  available  to  handle  the 
case  n=3.  A  collocation  analysis  for  n=3  would  be  feasible  but  would  be  rather 
cumbersome  and  would  not  necessarily  be  considered  reliable.  How  many  chord- 
wise  modes  would  one  have  to  use  with  spamvisc  waves  this  short?  From  this  point 
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of  view,  and  in  view  of  the  distinct  difference  in  the  character  of  the  chordwise 
distributions  in  (2.  2),  it  becomes  of  interest  to  take  a  look  also  at  the  pressure 
distributions.  This  we  will  do  in  Section  7. 

6.  Alternative  Elements 

An  alternative  to  the  set  wn(y)  *s  the  set  A*  which  was  defined  in  (2.  9).  If 

one  sets  II CL..Vl ,  =  w  ,  then  two  corresponding  elements  differ  by  the  effect 
4A'?+1  n,  n 

m  g 

of  Sj  according  U<  (2. 15)  resp.  by  the  contribution  of  c^>J(  according  to  (3. 1).  For 

the  element  A*,  the  set  €2^  is  immediately  given  by  (2.  7b).  The  downwash  w(y) 

which  belongs  to  A*  docs  not  fulfill  the  lip  condition ;  on  the  other  hand,  wc  saw  in 

Fig.  10  that  the  tip  singularity  becomes  less  and  less  visible  as  n  increases,  and  it 

is  illustrated  by  Fig.  12  that  the  two  elements  become  increasingly  more  equal 

over  an  increasing  range  of  the  wing  span. 

The  downwash  distributions  A*,  calculated  by  means  of  (3. 1)  and  (2. 14), 

are  shown  as  full  lines  in  Fig.  12.  All  A*  curves  turn  down  to  w(l)=-oo  but,  since 

the  downward  turn  occurs  closer  and  closer  to  the  tip  as  A*  increases,  a  finite 

limit  value  w(l-0)  is  built  rip  with  an  infinite  number  of  elements  A*.  Also  showr. 

are  (as  dashed  lines)  the  distributions  w  (y)  =  P-  ;  there  is  a  noticeable  difference 

n  411 

between  the  two  elements  when  n=0,  a  small  difference  when  n=l,  and  no  visible 
difference  (except  directly  at  the  tip;  not  showqlwhcn  n£2. 

The  set  A *  can  be  useful  because  of  its  analytical  simplicity  if  one  wants 
to  investigate,  for  example,  the  local  effect  of  a  disturbance  or  discontinuity  in 
w(y)  on  the  inner  part  of  the  wing  span.  In  such  a  problem  the  higher  w  plays  the 
important  role,  and  for  these  the  local  difference  between  the  two  sets  of  elements 
becomes  entirely  negligible.  Analytically  the  two  sets  differ,  of  course,  in  that 
the  A*  elements  do  not  exhibit  tho  complications  of  (2. 16);  however,  (2. 16)  describes 
the  tail  behavior  which  produces  the  proper  tip  singularity  in  the  n-elements.  The 
amplitudes  and  with  tf,A»  s  are  those  which  essentially  desribe  the  local 

pressure  distribution,  and  these  hardly  different  between  the  two  typos  of  elements. 


25 


7.  Pressure  Distribution 

The  pressure  distributions  in  Figs.  13  and  14  represent  the  cases  A*=l  and 
A*=3  but  would  not  look  noticeably  different  if  they  would  represent  the  cases  n=l 
and  n=3. 

Like  Fig.  2,  Figs.  13  and  14  show  relief  diagrams  of  p  over  the  right  half 
of  the  circular  wing.  The  left  half  of  the  wing  is  not  shown  in  Fig.  13;  it  is  shown 
in  Fig.  14  to  exhibit  the  wavy  distribution  of  the  wing  incidence.  In  both  figures, 
the  left  wing  half  is  stretched  into  a  rectangle  for  better  visibility  of  p;  furthermore, 
the  rectangle  is  cut  apart  at  y=0.  9.  The  finite  chord  at  the  wing  tip  which  is  pro¬ 
duced  by  the  stretching  operation  exhibits  the  kinked-linear  character  of  the 
chordwise  distribution  at  the  tip. 

Part  of  a  vertical  plane  through  the  stretched  leading  edge  is  shown  by  a 
line  of  constant  p  (the  scales  of  p  in  the  two  figures  are  arbitrary)  and  by  vertical 
lines  of  constant  y. 

The  most  striking  insight  of  technical  interest  from  Figs.  13  and  14  is 
that  the  waviness  of  w(y),  while  well  reflected  in  p.  in  act  ordance  with  (2.  9a), 
has  disappeared  almost  over  the  rear  part  of  the  wing  in  Fig.  13,  and  completely 
in  Fig.  14.  This  stabilizing  effect  is  the  cause  for  the  smallness  of  the  wave 
amplitudes  in  C  «  in  fig.  10.  However,  it  is  of  more  gciv.u,l  interest;  while  it 

^  i  n 

is  shown  here  on  idealized  models,  it  must  be  of  significance  for  the  influence  of 
any  disturbances  of  short  spanwise  extension. 

The  chordwise  distributions  of  p  in  Fig.  14  arc  not  kinked  (except,  of 
course,  at  the  tip  itself)  but  they  are  fairly  complicated  curves.  It  had  been  shown 
in  Ref.  8  that  for  the  planar  circular  wing  (n=0)  already  two  Birnbaum-Glauert 
chordwise  modes  yield  an  excellent  approximation.  For  the  case  n=3,  only  two 
modes  (i.  e. ,  approximation  by  a  parabola)  would  clearly  be  insufficient. 

To  the  degree  that  we  can  interchange  Xs:t=3  and  n=3,  the  difference  between 
the  complicated  distribution  Fig.  14  and  the  reference  distribution  Fig.  2  is 
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contained  in  the  column  C~  .  ,  of  Table  II.  This  remark  illustrates  how  the  bulk 

2h,  3v 

of  the  complete  set  describes  essentially  the  tip  singularity,  while  the  details 
of  p  over  the  inner  part  of  the  wing  are  fully  contained  in  the  much  more  rapidly 
converging  void  set  C ^  v» 

8.  Conclusion 

Since  the  nature  of  the  pressure  singularity  at  a  parabolic  wing  tip,  and  in 
particular  all  required  details  of  the  structure  of  exact  solutions  for  circular  wings, 
are  now  known,  it  has  become  possible  to  use  numerical  results  with  confidence. 
One  knows  about  their  convergence  behavior  and  can  match  with  the  asymptotic 
wing  tip  solution.  A  simple  procedure  for  calculating  numerical  solutions  for  the 
circular  wing  is  described.  It  is  so  far  specialized  to  symmetric  downwash  dis¬ 
tributions  without  camber  but  could  be  generalized  to  deal,  for  example,  with  an 
axisymmetric  fan-in-wing  configuration. 

Four  sample  solutions  have  been  calculated  to  very  high  accuracy  and  have 
been  tabulated.  The  solving  sets  converge  slowly;  their  bulk  describes  the  tip 
singularity,  typified  by  the  solution  for  the  planar  wing.  This  has  been  split  off; 
the  specific  solutions  which  remain  are  void  of  the  tip  singularity  and  converge 
rapidly.  The  sample  solutions  exhibit  details  of  the  mechanism  of  lifting  flow 
which  are  of  technical  interest. 

The  numerical  technique  which  is  described  in  this  paper  allows  one  to 
investigate,  for  example,  the  pressure  distribution  which  arises  from  a  local 
disturbance  or  discontinuity  in  the  downwash.  In  such  a  problem  approximation 
by  means  of  the  alternative  elements  of  Section  6  may  be  useful,  since  this 
approximation  becomes  increasingly  correct  with  higher  order  (i.  e. ,  s  large) 
effects  and  since  its  analytical  simplicity  may  allow  analytical  treatment  of  the 
given  problem. 
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Table  IL  Void  Solutions  (n 


1 , 2  and  3 ) 


(y)  =  (21y*  -  14y£"  +  1 )/ S  -  0.2795  066 

(y)  =  <42?y6  -  495y4  +  I35y2  -  5)/ 64  -  0.20268387 
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Fig.  3.  Chordwise  Position  of  Center  of  Pressure 
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Fig.  ?.  Elementary  Downwach 


Fig.  12.  Downwash  Distributions  A* 


APPENDIX 


INFINITE  SUMS 


All  the  elementary  sets  E^  from  which  the  solving  set  V  is  built  up 
obey  the  zero  sum  condition  (2.6);  the  leading  element  Eq  is  thus  given  by 


Ecr 


(A.l) 


Since  the  EK  are  reciprocal  progressions,  the  rate  of  convergence  of  the  sum  slows 
down  indefinitely  as  H.  increases,  and  numerical  term-by-term  summation  is  not 
usually  a  practical  proposition. 

Some  of  the  required  sums  are  tabulated,  see  e.g.  (2. 10)  for  the  rational 
sets  E^.  For  the  first  set  in  (2. 16)  we  have  from  6.3.  16  of  Ref.  6 


^j-TT  =  4(V(5/4)  +  ^]=  1.  399048526..  (A. 2) 


Certain  sums  involving  LK  are  derived  in  Ref.  5;  for  example 


(A.  3a) 


=  2. 103599581 


j>  ^  =  [31^(5)  -  14*(2)<f(3)]/4 

l  vv* 


=  1. 115624875.. 


(A.  3b) 


The  sum  (A.  3b)  occurs  in  Tables  I  and  II.  The  sum  (A,  3a)  occurs  in  the  leading 
term  of  C^+1,  see  (2. 16c). 

When  tabulations  are  not  available,  one  can  use  the  formula 


Vif  .*••• 


(A.  4) 


to  calculate  the  tail  sum.  Equation  (A.  4)  assumes  that  F,^  converges  smoothly 
and  that  the  function  E(x)  interpolates  E^  such  that  E(K)  =  E^.  Equation  (A;  4) 


W  -******' 


was  used  for  example  for  the  extrapolations  based  on  (4, 7). 

2 

For  the  term  with  log  k  in  (2. 16)  one  requires  the  sum 


*  0.  06505816 

nf 


which  was  determined  by  means  of  (A.  4). 


